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I. INTRODUCTION 



Recently the gauge/gravity dualityjl-3| has been widely applied to the study of con- 
densed matter physics. In particular, some critical phenomena in strongly coupled systems 
can be described by the dynamics of geometry as well as matter fields in a semi-classical re- 
gion (For recent reviews we refer to |4H7|). One remarkable example is building holographic 



mp- 
. In 



superconductor models based on the Abelian Higgs mechanism, through which an as 
totical anti-de Sitter black hole can break the U(l) gauge symmetry spontaneously [8 
these models, when the effective mass of a scalar field in the bulk is below the Breitenlohner- 
Freedman bound J9], a hairy solution to the scalar field can be obtained, characterized by its 
condensation around the horizon of the black hole . According to AdS/CFT correspondence, 
the creation of charged condensation strongly implies that a second order phase transition 
could occur in the dual CFT 



10 



13[| . Furthermore, this sort of holographic superconductors 



contain more interesting features 
the predictions of the conventiona 



14M25| . For instance, the energy gap is much larger than 



ductors as found in experiments 



26j 



BCS theory but quite similar to the high-T c supercon- 



In addition to various holographic s-wave superconductor models with scalar fields, one 
can also construct p-wave superconductor models with vector fields, implemented by intro- 
ducing a SU(2) Yang- Mills gauge field 27H32| . Correspondingly, the dual CFT has a global 
SU(2) symmetry and hence a conserved current J 1 . In the AdS black hole background, it 
is found that the U(l)s gauge symmetry (a subgroup of SU(2)) is possibly spontaneously 
broken such that the value of (J].) is not vanishing. This corresponds to a phase transi- 
tion between a non-superconducting state at high temperature and a superconducting state 
below the critical temperature. However, different from the s-wave model, here the order 
parameter is a current and the conductivity is anisotropic since the condensation of vector 
field breaks the rotational symmetry as well. As a result the conductivity has two indepen- 
dent components. One perpendicular to the direction of the condensation behaves like that 
in s-wave superconductors, while the other parallel to it performs a much different behavior. 
In particular, this component agrees well with the Drude model in the low frequency limit. 

Recently a gravity theory with nontrivial curvature- cubed terms in five- dimensional 



spacetime, usually called quasi-topological gravity, has been proposed in 



23-1221 . This the- 



ory can be viewed as a generalization of Gauss-Bonnet gravity. Besides the GB term in the 
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action, it also involves in higher-derivative corrections, thereby corresponding to CFTs with 



more couplings between operators as discussed in |38l . [39|. As a matter of fact, the holo- 
graphic study of the quasi-topological gravity has been carried out in many references 39M43| 
and its dual CFTs display much richer structures and novel features. Specially, in this the- 
ory two central charges relating to the conformal anomaly can be unequal, which brings in a 
non-zero but much lower bound of the ratio of shear viscosity to density entropy and hence 
violation of the Kovtun-Son-Starinet (KSS) bound 39(. In this paper, we intend to continue 
our previous investigation on holographic s-wave superconductors in 4l|, and construct a 
holographic p-wave model in the framework of quasi-topological gravity. We are specially 
concerned with the anisotropic behavior of the gauge fields and intend to compare our results 
with other p-wave models in Einstein's and Gauss-Bonnet gravity. We will also discuss the 
charge transport using a linear response theory, with a special interest in its behavior in the 
low frequency limit. 

We organize our paper as follows. In Sec. II, we present the holographic setup for a p-wave 
superconductor in quasi-topological gravity, then study the superconducting phase transition 
in the probe limit, focusing on the variation of the critical temperature with the coupling 
parameters. Sec. Ill contributes to the numerical evaluation of the anisotropic conductivity. 
Based on the data obtained we mainly discuss the following two issues. One is on the change 
of the ratio of the frequency gap to the critical temperature with the frequency, and the 
other is on the low frequency behavior of the conductivity. Discussions and conclusions are 
given in Sec. IV. 



II. QUASI-TOPOLOGICAL HOLOGRAPHIC P-WAVE SUPERCONDUCTORS 

We start with the five-dimensional quasi-topological gravity with an SU(2) Yang-Mills 
gauge field. The bulk action is given as 

S Mt = / JV^fij-^* + % + ^* + ^ ~ 4^^)] , (1) 

where G§ is the Newton constant in five-dimensional theory, and a, (3 and qym are Gauss- 
Bonnet coupling parameter, curvature-cubed interaction parameter and Yang-Mills coupling 
parameter, respectively. F l is the field strength of Yang-Mills gauge field with SU(2) gauge 
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symmetry and i is the internal index. Here A4 is the Gauss-Bonnet term 



A4 = R^ LVpa R II ' vpa — AR^R^ + R° 



(2) 



and Z5 is a curvature-cubed term with the form 

-£5 = R^u pa R P u a ^ R a ^ v + -yi^-R l ivp(jR iiV( " T R — 120R tlupcr R pup a R CTa 
+U4R pupa R flp R UCT + 128R/R/R P P - IWRJRSR + 11R 3 ). 



(3) 



It is worthy to point out that in contrast to higher order terms in Lovelock gravity [44|, the 
cubed terms above are not just topological but have contributions to equations of motion for 
bulk fields. Through this paper we will only take account of the probe limit of the theory. 
Namely, we will neglect the back reaction of the Yang-Mills field on the background metric 
in the large N c limit, where N c is the number of degrees of freedom per point in the dual 
free field theory 1 . 

In this limit stable AdS black hole solutions in five- dimensional spacetime have been 
found in 33| and they can be described as 

- 2 L 2 



ds 2 = _(-jV( r ) 2 /(r)dt 2 + dx 2 + dy 2 + dw 2 ) + 

L z r z j (r j 



where f(r) has three different solutions for different regions in parameter space 

a 



fi(r) = u + v - 



3/3' 



/2W 
Mr) 



--(u + v) 



•A ^ 

1 — (u — V) 
2 v ; 



a 

w 

a 



1 / \ ,v3 , . 

■-(u + v) — 1 — (u — v) 

2 V ; 2 v ; 3/3 



(4) 



(5) 
(6) 
(7) 



with 



u 



P 



(q+V7 z 

3/3 + a 2 
9/3 2 ' 



P 



3^/3 



2a 3 + 9a(3 + 27/3 2 (l - 
54/3 1 



(8) 



1 The large N c limit implies G5 ~ N c 2 — > 0, leading to a decoupling between the matter field with a finite 
Qym and the gravity in the action Sbuik- 
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L is the AdS radius and N(r) = N = \j \j fir) | r ^oo is the lapse function 2 . Now it is 
straightforward to obtain the Hawing temperature of these black holes, which is 

It will also be viewed as the temperature of the dual CFT on the boundary. 

Now we turn to construct the holographic p-wave superconductors. Firstly we need to 
solve the Yang-Mills equations in a fixed black hole background. Following the strategy 



presented in |8 



we take the ansatz as follows 

A a = Ay{dx») a = 4>(r)T 3 (dt) a + 4>(ry(dx) a , (10) 



where r % = a l /2i (i=l,2,3) with commutation relations [t\ t j ] = ie % i k r k are SU(2) gener- 
ators. In f fTUj) the nonvanishing vp(r) breaks the f (1)3 gauge symmetry generated by r 3 . 
We may interpret it as the p-wave superconducting phase transition from the side of CFT 
on the boundary, since in the dual field theory the global f7(l)3 symmetry is broken and 
superconducting charges can be created from the new vacuum which corresponds to the 
formation of the Cooper pairs. For convenience, we absorb the gauge couplings into the 
rescaling of the gauge fields 

4> = g Y AiL 2 (j), tp = gYAiL 2 ijj. (11) 

Moreover, we redefine the coordinate z = — = - such that the position of the horizon is 
fixed at z — 1, while the boundary is z — > 0. Then with the ansatz in Eq.f llOj) the equations 
for Yang-Mills field reduce to the following form 

r .l-^f, = , (12) 

z z l g 

r +(-+ V + -J^-^ = °> ( 13 ) 

g z i\ z g z z 4: 

where g = r ^ and the prime denotes a derivative with respect to z. Before solving these 
two equations we give the boundary conditions near the horizon and near the AdS boundary 
as follows: 



In order to get a normalized velocity of light on the boundary, the lapse function N usually should not be 
set to unit. This gauge is different from that in [3^] and some discrepancy of our results for GB gravity 
with those in 13211 c an be ascribed to this different gauge. Moreover, the argument that N(r) is a constant 
can be seen in |33f . 
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♦ The regularity condition at the horizon (z = 1) requires the gauge fields should be ex- 
panded as 

^ = ^ 0) +^ 2) (l-^) 2 + --- 

= <f>U(l- z ) + ... . (14) 

♦ Near the AdS boundary [z — > 0), the asymptotical behavior of fields are like 

^ = ^(o) + ^(V + --- 

= n + pz 2 + ■ ■ ■ . (15) 

In AdS/CFT dictionary, \i is the chemical potential on the boundary while pt = 2p and 
Pn — 4>h are understood as the total charge density and the charge density in the normal 
state, respectively. So the p-wave superconducting charge density is p s = p t — p n . A nonzero 
if)^ and nonzero i})^ correspond to a source and the expectation value of vacuum for the 
current operator J], that is dual to the gauge field A\ = ip respectively, so we have 

<4>=V> (2) . (16) 

For normalizable modes, the expectation value of vacuum can be obtained by setting ip^ = 
0. 

Before doing the numerical analysis, we present some remarks on the allowed range of 
the values of the coupling parameters a and (3. First of all, to obtain the stable black hole 
solutions without ghost modes or naked singularity, one finds that the allowed parameter 
range is constrained in the region as illustrated in the left plot of FIG Jt|33 ]. On the other 
hand, to ensure the positivity of central charges, energy flux and a well-defined causality for 
the dual CFT, a and /3 are further severely confined into a small region as dictated in the 



right plot of FIGJTJ which has originally been presented in 39(. Thus in our paper we will 
investigate the holographic superconductivity with parameter values restricted in this small 
region. However, in this region the coupling constant (3 is severely restricted in a narrow 
interval, roughly from (3 = —0.001 to /3 — 0.001. Graphically when one changes the value 
of f3 in this region, the plotting is probably not sensitive enough to illustrate the changes 
of the physical quantities with different values of the coupling constants. Thus we take two 
actions in our plotting. One is to enlarge a local region in figures to demonstrate the shifting 
tendency of the curves with the values of coupling parameters. Secondly, we also take some 
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FIG. 1: The allowed range of the values of coupling parameters a and j3: In the left plot, the yellow 
region is for the stable solution and the different regions divided by the color lines correspond to 
different stable black hole solutions given by (5) — (7) (see 3] for details); In the right plot, the 
blue region is the valid parameter range restricted by consistency of the dual CFT, and the red 
line corresponds to j3 = 0, i.e., the constraint of a in Gauss-Bonnet gravity. 

value for from the region in the left plot of FIG.l for comparison, for instance /3 = 0.1. In 
the probe limit we are allowed to do this since the back reaction of the perturbations is not 
taken into account and we may still obtain the stable phase in dual field theory, but it is 
warned that the curves obtained for these values are potentially instable for a system when 
the back reaction of the perturbations is considered, thus should not be trusted seriously. 

Now to explore the p-wave superconducting phase, we need find nonzero solutions for 
ip by numerical analysis. In our program, we set L — 1 and find the numerical solutions 
to the differential EOMS from the horizon to the AdS boundary, namely, from z — > 1 to 
z — y 0. In Figure El we illustrate the condensation of J\ as the function of the temperature 
with different values of coupling parameters a and /3. Note that different from the case 
of s-wave superconductor, here we plot the dimensionless quantity with respect to 

J- c 

^- since the conformal dimension of Jl is 3 rather than 1 for T c . FIG. 2 indicates that 

J-c x 

the order parameter has the behavior (J*) oc (1 — T/Tc) 1 / 2 near the critical temperature, 
and the value of \ x increases with both the Gauss-Bonnet coupling parameter a and the 
curvature-cubed coupling parameter /3, which is similar to the phenomenon obtained in the s- 
wave superconductor in quasi-topological gravity 4l| , as well as the p-wave superconductors 



in Gauss-Bonnet gravity 321 ]. 



In figure El we plot p s /pt v.s. by changing either of the coupling parameters a and 
p. The critical temperatures for different coupling parameters can be read off from the 



FIG. 2: The condensate as a function of the temperature with different values of coupling param- 
eters. In the left figure, the Gauss-Bonnet parameter a is fixed at 0.06. j3 is fixed at 0.001 in the 
right figure. In both cases the condensation tends to increase with the coupling parameters. 



u=0.06 P-0.001 




FIG. 3: The ratio of the superconducting charge density to the total charge density v.s. the 
temperature with different values of coupling parameters. 

intersects with the horizontal axis and their values are listed in TABLE I. We can see that 
the critical temperature decrease as either a or f3 increases, implying that with the increase of 
the curvature-cubed term or the Gauss-Bonnet term the occurrence of condensation should 



become harder. This tendency is the same as those in the previous papers[32|, |4l|. For 
explicitness we demonstrate the dependence of the critical temperature on either of the 

we know the conformal field dual to the 



coupling parameter in FIGIU Moreover, in 



quasi-topological gravity is characterized by central charges, c and a, and flux parameters, 
tl and £4. Explicitly, these parameters can be related to the coupling parameters in the 
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a 


p 


T c 


0.06 


-0.0008 


0.1924/j 1 / 3 


0.06 


0.000001 


0.1923/9 1 / 3 


0.06 


0.001 


0.1921/) 1 / 3 


0.06 


0.1 


0.1667/9 1 / 3 


-0.2 


0.001 


0.2218/j 1 / 3 


-0.1 


0.001 


0.2120/9 1 / 3 


0.0001 


0.001 


0.200V/3 


0.1 


0.001 


0.1858/9 1 / 3 



TABLE I: The change of the critical temperature with the coupling parameter a and /3. The critical 
temperatures corresponding to the parameter values in FIG J3] are listed in two tables respectively. 




FIG. 4: The relations between T c and the coupling parameters a and /3, respectively. In the left 
figure the line is for a = 0.06, while in the right figure the line is for /3 = 0.001. 



bulk as follows 



c-a - 3/3/qq) 

P2 ' \ ') 



1 - 2^/00 - 3/3f ( 



oc 



_ 24/00(0-87/00/3) 

1 - 2«/oo - 3/3/2 



' 00 

C2 



3780/4/3 

" l-2a/oo-3/3/^' (1J) 

where /oo satisfies 1 — /oo + cc/^ + /3/^, = 0. Among these parameters any two of them is 
enough to describe the correlations in the dual CFT. For example, we can calculate 5 and 
£4 as the free parameters and fit how T c changes with the changing of 5 and tA respectively. 



9 



However, as discussed in [4l|, a clear rule from the CFT side is still missing, thus we will 
not show the fitting plot of the critical temperature versus 5 or t4 in our current paper. 

In the next section, we turn to investigate the conductivity and find its new characters 
comparing with the s-wave superconductors. 

III. CONDUCTIVITY 

In this section, we will study the charge transport and linear response of the boundary 
system. In the linear response theory, a central quantity is the retarded Green function. 
According to the AdS/CFT dictionary, if we want to know the retarded Green function of 
the U(l) current, we just need to study the propagation of the linear perturbation of the 
gauge field in the bulk in the probe limit. So, in the following, we are interested in the 
linear response of the r 3 component of the Yang-Mills gauge field. As discussed in jsj], we 
can consider an alternating current (AC) on the boundary by introducing a time-dependent 
perturbation for the gauge field 

A->A + 5A, (20) 



where 

5 A = e~ iu>t 



a\{r)T l + a 2 t (r)T 2 ^jdt + a 3 x (r)r 3 dx + a 3 y (r)r 3 dy . (21) 



Though the condensation of ip breaks the rotational 5*0(3) symmetry associated with x 
direction, the system still has an SO (2) symmetry in y — w plane. Hereafter, we neglect the 
perturbation along to- axis and only consider the electrical conductivity a xx and o yy . 

1. Oyy component 

Substituting the ansatz in (1211) into the Yang-Mills equations, we obtain the equation of 



motion for a 3 



s" , ,9' , In 3' , ( w L w \ 3 



where the prime is relative to z — 1/r. In general, different modes of the perturbations are 
mixed, however, from the equation above we notice that a 3 is actually decoupled from other 
components. In addition, the equation is similar to that for s-wave superconductors 
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4l|, 
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so we can obtain the conductivity a yy in a parallel way. Specifically, we choose the ingoing 
wave condition for a 3 , near the horizon, then we have 

oJ = (1 - z)&[l + a f\l -z) + af\l -zf + ■■■}. (23) 

Moreover, the behavior of a 3 in the asymptotical AdS boundary (z — > 0) is 

afViV 4 



a 



I = af ) + af V + "» | (logA/,), 2 , (24) 



where af , af and A are integration constants. The last term in ( l24"j) will lead to a diver- 
gence when one calculates the Green function, however, such a logarithmic divergence can 



be canceled with a boundary counterterm in the renormalization procedure 45[ . Thanks to 



the standard AdS/CFT dictionary, the conductivity can be expressed through the retarded 



Green function as follows [461] 



a(u) = —G r (oj) 

lOJ 



— lim Ng(r)ra 3 a 3 ' (25) 



k=0 IU) r-»oo y v 



Thus we find the conductivity a yy is 



2af } 



a 



VV ' ■ AT3 3 (°) 

iu)N 6 a y 



iuNlnN + . (26) 



Given the boundary conditions ( |23l) and (l2~4l) . we numerically solve equation ( |22l) to obtain 
af •* and af \ As a result, the dependent relation between a yy and u; for different coupling 
parameters /3 and a are shown in FIG JS] and FIG® respectively. From these two figures, it 
is clear that when the frequency u is large enough, the real part of a yy always increases while 
the imaginary part becomes linear due to the dominant term uNlnN + ^-. In addition, from 
the imaginary part of conductivity one finds that the ratio ou g /T r increases with either a or 



(5. This behavior is much similar to that of a in the s-wave case|41|, while it is interesting 
to note that for p-wave superconductor the values of the real part of conductivity increase 
with the coupling parameters. Since this quantity corresponds to the imaginary part of the 
Green function which characterizes the dissipation of the charge transport, this phenomenon 
is analogous to what happens in holographic hydrodynamics as described in figure 2 in 



3 As pointed out in 47[, the term —iuNlnN should not be neglected. Though we take a different gauge 
here, the asymptotic analysis near the boundary is similar and a direct calculation shows that the result 
is the same 
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|5=-O.O0O8 
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J5=0.001 



t 0.04 - ~ ■ 
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5=0.001 
W).1 
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FIG. 5: The conductivity a yy for the p-wave superconductors with a fixed Gauss-Bonnet parameter 
a = 0.06. The left figure is for the real part of a yy while the right one is for the imaginary part of 

<Jyy. 

P=0.001 P=0.001 



I 


u=-0.2 * 
-a=-0.1 , ' 




- a=0.0001 * " 




~ 0=0.1 
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FIG. 6: The conductivity a yy for the p-wave superconductors with a fixed the curvature-cubed 
parameter f3 = 0.001. Similarly, the left figure is for the real part of a yy while the right one is for 
the imaginary part of a yy . 

2. a xx component 

Now we intend to study the component a xx , which can be obtained by solving the equation 
of motion for a x in the bulk. In contrast to dy, since the component a x couples with the 
fields a\ and in the linearized Yang-Mills equations, besides three second-order coupling 
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equations of motion 



s" A d\ 3' o iuaf + a\4> , 

a ? _ laf + Ef* = o, 

z z 2 g 



af - "-a* - ^a 2 - = 0, (27) 

z z A g z z g 



L 2 ip 2 2 iL 2 ua x ip 
YT a t IT, 



we need to solve another two first-order equations together 

— iua\' — (pa 2 ' + a 2 (p' = 0, 
cj,a]' - iuaf + L 2 N 2 g^z 2 - a t V - L 2 a\N 2 gi,'z 2 = 0. (28) 

Again using the ingoing wave condition near the horizon, we have the following asymptotical 
behavior of a^, a} and a 2 : 

► Near the horizon (z — > 1) 

4 = {l-z)&[l + Af\l -z) + Af\l -z) 2 + ■■■}, (29) 

a] = (1 - z)^[Af\l - z) 2 + Af\l -zf + ■■■}, (30) 

a\= (1 - ,)w[ A 2(1) (l - 2 ) + A 2(2) (l - ,) 2 + • • •]. (31) 

► Near the boundary of the AdS bulk {z — > 0) 

al = af ] + af\ 2 + ^ 
a] = af ] + af ] z 2 + --- , (33) 

a 2 = G 2(°) +G 2( 2 ) z 2 + ... _ (34) 

All the coefficients can be determined numerically. However, these quantities are gauge 
dependent since we have not done any gauge fixing to these components. As discussed in [8j , 
in order to obtain the conductivity a xx we need define a gauge invariant quantity as follows 

- 3 iuL 2 a 2 + (pa\ 

al = a x+ 02-a;2jL4 (35) 
Then we have the asymptotical behavior of a x near the AdS boundary 

o(o) ,(2) o al i0) u 2 NHog(A/z)z 2 iuL 2 a 2i0) + /ia, l(0) , 2) 2 
a3 = a 3 x { ' + af ' z 2 + -2 „ + j TTJ^^ z 

^ af > + ^V + 4'V^MW (36) 
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FIG. 7: The conductivity a xx for the p-wave superconductors with a fixed a. The left figure is for 
the real part of a xx while the right one is for the imaginary part of a xx . 



- o=-0.2 

- a=-0.1 

- u=0.0001 



"Hi' 
iii i - T - - 

til N 



20 25 30 35 



FIG. 8: The conductivity a xx for the p-wave superconductors with a fixed curvature-cubed param- 
eter j3 = 0.001. Similarly, the left figure is for the real part of a xx while the right one is for the 
imaginary part of a xx . 

where we have defined 

ln\ , * , T 2„2(0) j_ ,,„i(°) 

af = af } + ^ ? + ^ (37) 



As a result, the conductivity a xx has the form 

1 2in 3 ^ 

— Ct [UJ k =0— TTTyj — 



a ra = —G R (oo) | fc=0 = - * - iwMniV + ^iwiV. (38) 



In order to understand the behavior of a xx , we simultaneously solve the equation groups 
(P|b (fl3"|h p7|) as well as (j25). The numerical results of are shown in FIGdand FIGJHl 
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B -7 -6 -5 _4_3_2-l 1 

log ]0 (iu/p A 1/3) 



-8 -7 -6 -5 -4-3-2-1 1 

log 10 (w/p A 1/3) 




-8 -7 -6 -5 -4 -3 -2-1 1 -8 -7 -6 -5 -4 -3 -2-1 1 

log 10 (M/p n 1/3) log 10 (w/p n 1/3) 



FIG. 9: The relations between logio[Re(a /T c )] and logio[u/ // 1//3 )] for fixed a = 0.06. The blue 
and black lines are for a yy and a xx respectively, while the red ones are the fitting lines of a xx when 
to is low enough. 



a 


P 


do/T 


rT 


0.06 


-0.0008 


1.0213 x 10 9 


1.8304 x 10 6 


0.06 


0.000001 


7.5091 x 10 8 


1.3064 x 10 6 


0.06 


0.001 


6.4602 x 10 8 


1.1072 x 10 6 


0.06 


0.1 


5.6827 x 10 8 


6.38 x 10 5 


-0.2 


0.001 


3.735 x 10 9 


8.804 x 10 6 


-0.1 


0.001 


2.5534 x 10 9 


5.588 x 10 6 


0.0001 


0.001 


1.5959 x 10 9 


3.2288 x 10 6 


0.1 


0.001 


8.7646 x 10 s 


1.6036 x 10 6 



TABLE II: The evaluation of ctq and r and their dependence on the values of coupling parameters. 
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-9 -8 -7 -6 -5 -4 -3 -2-1 1 -9 -8 -7 -6 -5 -4 -3 -2-1 1 

log 10 (oi/p A 1/3) log 1f) (o)/p A 1/3) 



FIG. 10: The relations between log 10 [Re(a/T c )} and log w [u / ' p^] for fixed = 0.001. The blue 
and black lines are for a yy and a xx respectively, while the red ones are the fitting lines of a xx when 
uj is low enough. 

In these two figures, we find that the point of phase transition will also change when we 
choose different couplings, but the behavior of the component a xx is very different from that 
of component a yy . Comparing with a yy , the real part of a xx ascends much more slowly. In 
the high frequency limit, the imaginary part has a suppression effect with the increase of the 
coupling parameters. On the other hand, in the low frequency limit, a xx behaves quite like 
that in Drude model, and this analogy was firstly observed by Gubser in Drude model 
is conventionally used to describe the classical electron system, in which the real part of the 
conductivity is featured by the relation 



Re(a) = — ^-j, (39) 

where o"o = ne 2 r/m is the DC conductivity, n, e, m, r are the electron density, charge, mass 
and the relaxation time respectively. Based on the Drude relation, we can fit our numerical 
results to evaluate o"o and r in FIG. [9] and FIG fTUl All the fittings are done under the 
condition that the condensation is stable, namely, the value of T/p 1 / 3 is small enough. We 

1(3 



show the fitting results of o"o and r and their dependence on couplings parameters a and 
in TABLE II. From the table, we find that both of o"o and r depend on the coupling 
parameters. 



IV. DISCUSSIONS AND CONCLUSIONS 



In this paper we have constructed a p-wave holographic superconductor model in quasi- 
topological gravity in the probe limit. Firstly, we find that the superconducting condensation 
becomes harder with the increase of the GB coupling and curvature-cubed coupling. This 
is very similar to the case of the s-wave model. Secondly, both anisotropic conductivities 
a xx and a yy are studied. The numerical results indicate that the behavior of a yy at the low 
temperature is quite similar to that in the s-wave model. More precisely, with the increase 
of coupling parameters a or (3 the ratio oj 9 /T c m 8 becomes unstable and increases as well. 
However, the conductivity a xx behaves differently. Its real part grows more slowly with the 
frequency, while its imaginary part contains a spike. For both components, the imaginary 
part is always suppressed by the increasing of couplings in the large frequency limit. While 
in the low frequency limit, our data of a xx fits the Drude model very well, but the values of 
DC conductivity as well as the relaxation time depend on the coupling parameters. 

It is worth pointing that due to the presence of higher curvature corrections in quasi- 
topological gravity, besides the ghost modes and naked singularity one should also be cau 



jility as discussed in [48] 
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491 1 . Those potential 



tious of other potential instabilities, such as the dynamical insta 
for Lovelock gravity and the plasma instability as discussed in 

instabilities may further restrict the valid values of the coupling parameter s [39j, |48|, |49j. We 
leave this open issue for further investigation in future. 

In the end of this paper we remark that it should be very worthy to investigate p-wave 
superconductors when the back reactions are taken into account in our model. It is expected 
that both the condensation and the charge transport would be corrected by the effects of back 
reactions 
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30 



32 
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501 ] . Moreover, inspired by recent progress on the holographic non- 



fermion liquid and strange metals 



51 



U456I]. it might be possible to explore fermion system 
with a finite charge density in the framework of the quasi-topological gravity. In these 
systems, the fermion part near the horizon has a loop contribution to the total two point 
correlator of the current, however, such an 0(N°) contribution dominates the dissipation 



17 



of the transport [571] . Since in quasi-topological gravity the bulk of spacetime may exhibit 
a richer structure of geometry due to the higher order couplings, we propose that the dual 
CFTs would show different behavior at low energy limit. This is expected to be done in 
future. 
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